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Abstract 

The resonance modes and the related effects to the transmission of elastic waves in a two dimen- 
sional phononic crystal formed by periodic arrangements of a two blocks unit cell in one direction 
are studied. The unit cell consists of two asymmetric elliptic cylinders coated with silicon rubber 
and embedded in a rigid matrix. The modes are obtained by the semi-analytic method in the least 
square collocation scheme and confirmed by the finite element method simulations. Two resonance 
modes, corresponding to the vibration of the cylinder along the long and short axes, give rise to 
resonance reflections of elastic waves. One mode in between the two modes, related to the oppo- 
site vibration of the two cylinders in the unit cell in the direction along the layer, results in the 
total transmission of elastic waves due to zero effective mass density at the frequency. The reso- 
nance frequency of this new mode changes continuously with the orientation angle of the elliptic 
resonator. 

PACS numbers: 43.40. +s, 46.40. Cd, 63.20.-e 
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I. INTRODUCTION 



The propagation of elastic or acoustic waves in periodic heterogeneous materials has re- 
ceived renewed attention in the last several years^ ^^^^^ . A new class of composite 
materials with periodic structures in the elastic properties has been proposed to realize the 
sound attenuation at selective frequency bands in the audible range, from several hundred 
hertz to few kilohertz 1 . For suitable periodicity and elastic contrasts between the compo- 
nents, low-frequency gaps may appear in these so-called phononic crystals (PCs) in the form 
of a thin layer, which can lead to promising applications in sound insulation^. The periodic 
structures can also be applied to the design of acoustic wave guides or filters^. 

It has been demonstrated the improvements on the sound insulation of the PCs in com- 
parison with the homogenous materials whose ability of sound insulation is governed by the 
mass law&£. Specific efforts were made to the identify and application of the PCs' character- 
istics in the low-frequency range. The so-called first resonance frequency, the local vibration 
of the constitution elements of PCs, which is usually a coated sphere or other simple shapes, 
was thus extensively examined^ 1 ^. 

Theoretical methods have been developed in order to understand the properties of elastic 
response of PCs and predict new phenomena. A conceptually simple method named plane 
wave expansion method (PWE), which treats the wave equation in the Fourier spaced, has 
obtained a vast of important results of the PCs. However, when used in problems with large 
elastic constant contrasts, it suffers from slow convergence and heavy computation efforts. 
On the contrary, the multiple scattering method^, looks a little bit complicated, overcomes 
these difficulties. The problem with it is that the method, though very efficient in simple 
elementary geometries such as spheres in three dimensional and cylinders with circular cross 
section in two dimensional problems, is not efficient to deal with nonelementary resonator's 
geometries. In this case, it is practical to resort to some computational techniques, such 
as Multiple Multipole method (MMP)^ 1 ^, finite difference algorithms^, lumped mass 
method (LM)£, and variational method (VM)£, etc. 

In addition to the above mentioned successful methods, mechanical models have been 
developed for cylindrical elastic resonators (2D PCs). The models are based on a point 
mass connected as a pendulum or by springs to a rigid or elastic matrix 2 ' 8 and are able to 
predict the dependence of the resonator on the material constants. These models are useful 
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to gain insight into the underlying physics. Moreover, the line profile of the transmittance 
spectra for the first resonance peak was analyzed by C. Goffaux— along this line. Their 
mechanical model, taking into account the interaction between propagating waves in the 
matrix material and local resonance effects, evidenced the equivalence of resonance scattering 
in PCs to Fano's interference phenomena. Quite recently, complete modal analysis has been 
performed for single cell with Finite Element Method (FEM). The resonance frequencies of 
modal analysis revealed by FEM are compared with the simulated wave attenuation peaks. 
A small phase shift was observed and interpreted as Fano type interferences^. However, 
since the methods are based on PC whose repeating unit cell contains only one resonator 
connected to a matrix, thus can not treat properties that associated with local couplings 
between two nearest resonators, symmetric or asymmetric. 

In this paper, we use a different approach to study a 2D ternary PC with asymmetric 
resonators. We investigate resonance modes resulting from PC with two connected res- 
onators as the unit cell, using analytic formulation combined with Least Square Collocation 
Method (LSCM)^^, whose more advanced version is known as multiple multipole method. 
Then the model is simulated by FEM to confirm the existence of these resonance modes. 
With asymmetric arrangement of the two resonators, we identified a new resonance mode 
attributed to the coupling effects of the two resonators in a unit cell. Wave transmission 
spectrum is also calculated by FEM. The line profile exhibited in transmission spectrum 
corresponding to this new resonance mode is also discussed. 

II. THE SEMI- ANALYTICAL TREATMENTS AND RESULTS 

The resonance modes of a spherically symmetric elastic resonator and its cylindrically 
symmetric counterpart in two dimension (2D) have been derived analytically under the 
assumption of rigid core and matrix— in 2005. In 2006, Hirsekorn et extended the treat- 
ment by removing the rigid assumption and employing the physical material properties. 
And the analytic model were confirmed by the simulation results. Here, we derive an ana- 
lytical formulation describing an asymmetric 2D model, combined with LSCM for boundary 
conditions. 

The model system we studied is the PC with the unit cell as an asymmetric arrangement 
of two elliptic hard cylinders coated with soft silicone rubber and embedded in a rigid matrix, 
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FIG. 1: The unit cell of the PC in this study, constructed by two blocks of the ternary material. 
The two blocks are placed with mirror symmetry about the middle plane. The side length of each 
block is d = 16mm, the outer radius of the coating is c = 7.5mm, short and long axis of the 
elliptical cylinder are a = 3mm and b = 5mm, respectively. 9 is the vibration angle, gives the 
polarization of the vibration with respect to the x axis, and 6q is the orientation angle, specifies 
the angle between the long axis of the upper elliptic cylinder and the x axis. 

the cross section of the unit cell is shown schematically in Fig. [TJ The region 1 is the rigid 
epoxy matrix, region 2 are the two soft silicone rubber coatings with their outer shape as 
circles, and region 3 are two hard cylinders of elliptic shape. The geometric and physical 
parameters are given in the figure caption. 

Suppose a long wave-length elastic wave is traveling along the x — y plane, or, r — <fi plane, 
perpendicular to the axis of the cylinder. The cylinder is vibrating around the cylinder axis. 
Consider the lower cell of the pair shown in Fig. [XJ, and suppose that the vibration angle is 
9 with respect to the positive x direction (the arrow direction in Fig. [I]), the cylinder will 
translate as a rigid body in both x and y direction, 

— m 3 uj 2 U cos 9 = J J[n r cos (f)T rr + (n^ cos <p — n r sin </>)r r — sin (pr^rdldcj) (1) 

and 

— m 3 uj 2 U sin 9 = / / [n r sin <pT rr + (n<f, sin <fi + n r cos (f))^ + cos (^T^rdldcf), (2) 
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respectively. And its rotation described by 

- I 3 uj 2 (j)o = r{n r T r(j> + n^T^rdldcj), (3) 



where U is the displacement of the cylinder, O is the amplitude of the rotation, m 3 is the 
mass of the cylinder, expressed as p 3 irabl, J 3 is the inertial of rotation with respect to the 
cylinder axis, given by p 3 lnab(a 2 + b 2 )/4:, with p 3 being the mass density of the cylinder, a 
and b are the short and long semi-axis of the elliptical cross section of the cylinder, and I 
is its length. Here, r rr and r r $ are components of stress tensor in medium 2 in a cylindrical 
coordinate system, and n r and n§ are the r and components of the outer normal unit 
vector on the cylinder surface, which are functions of (r, 0). The surface integration on the 
right hand side of (CQ) and (j2j) are the forces exerted on the cylinder from stresses in the 
coating in x and y direction, and the right hand side of ([3]) is the toque on the cylinder. 
The displacement field in medium 2 is described by the elastic wave equation 2 ^ 

(A 2 + 2/x 2 ) V(V ■ u) - /x 2 V x V x u + p 2 oo 2 u = 0, (4) 

where A 2 and p, 2 are the Lame coefficients , and p 2 is the mass density of medium 2. Mirror 
symmetry about the x — y plane indicates that the displacement in the lower block is just 
the mirror reflection of the upper block. This means that the y direction displacement of the 
reflection plane is zero. Thus we only need to obtain the displacement fields in one block, 
which we choose the upper block to solve. The equation of the displacement field in medium 
2 may be expressed in terms of potential functions as^ 

u = V$ + V x (Ve z ), (5) 

where $ and \l/ are scalar potential functions, satisfying the following scalar wave equations: 

V 2 $ + a 2 $ = 0, (6) 
V 2 ^ + /? 2 ^ = 0, (7) 



here a = tuy'p 2 /(A 2 + 2p 2 ) and (3 = cuyp 2 //i 2 , denoting the wave numbers of the shears- 
wave and compressional(P)-wave, respectively. The solutions for $ and \l/ may be written 

as 

oo 

$ = Y,l A k-Jk(ar) + B k N k (ar)](C k cosk<f) + D k smk<P), (8) 

k=0 
oo 

* = Y,l E MPr) + F k N k (f3r)}(G k smk<P + H k cosk<P), (9) 

fc=0 



where Jk{x) is the Bessel function and Nk(x) the Neumann function. The angle part in (JHJ) 
and ([9]) can be expressed in linear combinations of sin0 and cos0, while the single- valued 
requirement indicates that k is an integer. For practical reasons, the expansion of $ and ^ 
have to be truncated after K terms, where index fee {0, 1, • ■ ■ , K}. The displacement field 
is continuous across the boundaries between medium 1 and 2 (<9r 12 ) 

u r \ r=c = u cos 0, (10) 
w</,| r = c = —u sin 0, (11) 

and between medium 2 and 3 (<9r 23 ) 

u r\dr 23 — t/cos#cos0 + U sin 6 1 sin — r(l — cos0 o ), (12) 
M^,|ar 2 3 — ~ U cos 6* sin + [/ sin # cos + r sin (f> , (13) 

where u is the displacement of the epoxy. 

In a local cylindrical coordinate system (r, 0, 2) , the strains due to a displacement u are 
expressed as^ 

6 rr = ( 14 ) 

e ^ = 90 + aT' (15) 

1 fldu r du<t> uA 
^ = ^ = 2{r^ + ^-V)- (16) 

All other components are zero since they involve the u z component or cross derivatives with 
respect to z. The stresses are linearly related to the strains by^ 

T pq = X5 pq J2en + 2fie pq , P,ge{r,0}. (17) 

i 

The displacement field is determined when the coefficients A^-Hk are known, which then 
can be obtained from the boundary conditions, (fl2"|) and ([TBI . Since the coordinate system 
used is not coincident with the boundaries of the core cylinders, the wave equations ([6]) and 
(J7J) are not separable on the boundaries. In principle, the boundary condition should be 
satisfied on every point of the boundaries, in practice, the expansions of (jSj) and are 
truncated at Kt\i term so that only a finite number of independent equations are needed 
to solve for the unknowns. The simplest way to solve the problem is to choose as many 
representative matching points on the boundaries as the number of unknown coefficients and 



solve for the coefficients from boundary conditions on the matching points. However, the 
solution obtained this way is not necessary to satisfy the boundary conditions at other points 
on the boundaries, and worst of all, the equations resulted in this way may not be solvable. 
For example, the resulted equations may be singular, or even may not independent. This 
problem is solved by the so called LSCM. In this method we choose the number of matching 
points much larger than the number of unknowns and solve the over determined equations 
in the least square sense. The boundary conditions on the matching points are not exactly 
satisfied, but the overall deviations from the exact boundary conditions will be minimal and 
the solution is "smoother" in between matching points 2 --. The displacement U and u are 
predefined symbol constants, thus in our specific problem we have 8(K + 1) independent 
unknowns if the truncation of the expansion is K. On the boundary between region 1 and 
2, we choose L matching points and on the boundary between 2 and 3, M matching points 
are chosen, hence 2(L + M) equations. The inequation 2(L + M) > 8(K + 1) should be 
satisfied. 

With the obtained coefficients, the r rr and t t ^ can be obtained, and then the right-hand 
side of equations of motion ([I])-®. Finally, the U cos9, Usin9, and O can be computed 
from the equations of motion, in terms of u. 

We calculated the displacements and other physical properties of the model for #0 = 45°, 
the truncation K = 30, and the matching points L = 180, M = 180, and other parameters 
are given in the caption of Fig. [TJ Now we discuss our results in detail. A resonance 
rotational mode is found at / = 258Hz, denoted as mode 1. In this mode, the inclusion 
cylinders rotationally vibrate within the coating, while the epoxy matrix remains nearly 
stationary. This resonance mode could be described by a model with the shear deformation of 
coating^. However, according to the rule proposed by Wang et aU&, it makes no contribution 
to wave insulation. The reason for this is that the forces to the hosting structure contributed 
by the oscillator is zero in this mode. Fig. [2] shows the displacement of the core cylinders. 
It is observed that there are two resonances around / = 313Hz and / = 363Hz, denoted as 
mode 2 and mode 4, at these frequencies the amplitude of the displacements of the cylinders 
in both directions are very large and experience a very sharp (~ 71") change of phase. By 
solving for the vibration angle 9, it is found that the lower frequency resonance vibration 
is along the long semi-axis of the elliptical cylinders and the higher frequency resonance 
vibration is along the short semi-axis of the elliptical cylinders. 
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FIG. 2: The frequency dependence of displacement for rubber-coated elliptical cylinders in both 
x and y directions, denoted by solid line for U x = UcosO, and dashed line for U y = U sin 0, 
respectively. 

The appearance of two different resonance frequencies is not only a result of the variation 
of the coating thickness, it also relates to the contact area between medium 2 and 3 in 
vibration direction. Along the short semi-axis the effective thickness is larger than the 
one along the long semi-axis, however, the contact area is much bigger for short semi-axis 
vibration. These give rise to a higher resonance frequency in mode 4 than mode 2. A detailed 
description about the two modes without coupling effects can be found in [3}. Between the 
two resonance frequencies, resulted from the splitting of the degenerate modes in symmetric 
resonator case, the phases of the vibration in x and y direction are reversed. This phase 
reversion brings about an interesting mode, denoted as mode 3, which demonstrates some 
notable physical properties. 

To find out the implication of the phase reversion, we calculated the frequency dependance 
of force on epoxy matrix by the coated cylinder, 



F x = J J (n r cos 4>T rr + (n<£ cos — n r sin 0)r r?i — n§ sin (pr^cdldcf) (18) 
F y = (n r sin 0r rr + (n^ sin <fi + n r cos 0)r r< ^ + cos (jjr^cdldcj), (19) 



plotted in Fig. [3j It can be seen from Fig. [3]^a) that between the two resonance frequencies 
the force in x direction is gradually changed from positive to negative, while the force in y 
direction is always positive. The cross point frequency, denoted by f c , where the force in x 
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FIG. 3: The frequency dependence of (a) force acting on epoxy matrix by coated resonators in 
x and y directions; (b) displacement of cylinders in medium 2, which is an enlargement of the 
resonance part in FigfSJ The line with solid squares and line with solid trigons represent x and y 
direction, respectively. Both the force and displacement are normalized to the displacement u in 
medium 1. 

direction is zero, is close and a little higher than 332Hz. On this frequency, the amplitude 
of the displacement field in y direction is a minimum, as is shown in Fig. E](b). Here the 
effective mass density (EMD) for the coated resonator is a tensor and its a;-component is 
denned by F£ 3 = — pf 3 V23^ 2 M, where V23 = vrc 2 / is the volume of the coated cylinder. At this 
frequency, the EMD for the block defined 19 as p x e = <pipi + (02 + $3) P235 wno normalized by 
Pe — <PiPi + 02P2 + 03P3, is 0.06. The zerolike EMD for the block can give rise to an irregular 
wave propagation behavior in the x direction, which will be discussed in detail below. This 
irregularity can not be simulated by single resonator in one unit cell. The nearest neighbor 
coupling effect of resonators plays an important role in this resonance mode. 

Both P-wave and .S-wave can exist in solids. Generally, the energy will be trapped by 
excited resonators, then gradually passed to the surrounding material in the process if the 
matrix's absorption coefficient is not zero. In mode 2 and mode 4, resonance appear in x 
direction, the EMD tends to be infinity, thus elastic wave can expected to be totally reflected 
by non-dissipative material. Conversely, in mode 3, the EMD tends to zero in x direction, 
so significant elastic wave transmission is expected. 
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III. FINITE ELEMENT METHOD SIMULATION 



To validates the results, the finite element commercial software, COMSOL Multiphysics 
is used to simulate the vector displacement field in a unit cell with a pair of long coated 
elliptic cylinders immersed in epoxy, whose cross section is shown in Fig. [TJ Periodic 
boundary condition is used in the y direction. The material component properties and 
their corresponding dimensions are consistent with calculations in Sec. [Til Instead of rigid 
assumption for the matrix and inclusion, some typical material parameters are assigned to 
the matrix and inclusion cylinders: the Young's modulus for each components are, E\ = 
4.35GPa, E 2 = 117500Pa, E 3 = 16GPa, the Poisson's ratios are, v x = 0.368, v 2 = 0.469, 
u 3 = 0.44, and the densities are, p\ = 1180kg/m 3 , p 2 = 1300kg/m 3 , p 3 = 1134:0kg/m 3 . 

A. Resonance Modes 

Triangular paver elements were used to construct finite element mesh of a unit cell with 
elliptic inclusions with 11173 nodes and 21984 elements. The mesh size is usually set ac- 
cording to the shortest wave length expected during the events. For example, if N ew is 
the number of elements per wavelength required for accurate modeling, C m i n is the slowest 
wave speed, and f max the largest frequency experienced in a frequency sweep, then the mesh 
size is determined by A min = C m i n /(N ew f max ), (e.g. N ew = 6 for quadratic element shape 
functions and N ew = 10 for linear element shape functions). When rubber like materials are 
employed, the shear wave speed is typically the smallest and governs the mesh size needed 
in the solid. In the frequency range considered in this study, the wave-length in the matrix 
material is much longer than the side length of the unit cell. 

The upper and down side were set as periodic boundaries to represent the infinite extent 
in y direction, therefore the displacement of epoxy matrix in this direction is zero. The left 
and right edges were set as free sides receiving and transmitting traveling waves. 

Fig. H] shows the resulting simulated modes. For application reasons, higher rotational 
modes of the coating with an increasing number of integer nodes are not discussed here. 
Each one of the four modes related to the motion of the resonator, so all of the resonance 
frequencies of these modes are expected to be in the interesting low frequency regime. In 
mode 1, the core is rotating in the inclusion, inducing no net forces acting on the matrix, 
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FIG. 4: Vibration modes in 2D ternary PC with asymmetry resonators at four frequencies: 
f=258Hz, 313Hz, 332Hz and 363Hz. The direction and length of the arrows represent the direction 
and amplitude of displacement vectors, respectively. 
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thus make no contribution to the wave insulation. Mode 2 and mode 4 are formed by 
the substitution of the circular shaped resonators by elliptic ones. Thus the degenerate 
resonance frequency in circular case is split into two modes. The resonance frequency in 
mode 2 is / = 313Hz, lower than the circular one (a = b = 5mm), which is / = 359Hz for 
the parameters used here, according to Liu et alM. The resonance frequency of mode 4 is 
/ = 363Hz, which is higher than the circular counterpart. 

Mode 3 is the result of the interaction between two blocks in the unit cell, the resonance 
frequency is f — 332Hz. The resonators are placed at geometry centers of the two blocks, 
while the orientation angles of the two resonators in one unit cell possess mirror symmetry 
with respect to the middle plane (see Fig. [TJ. The conjugate orientation angles for upper 
block and lower blocks in unit are 45° and —45°, respectively. The coupled two blocks give 
rise to a vibration mode, in which two resonators move in opposite directions along layer 
extension. The orientation angle can be set in the range 9$ G (0°, 90°) for the upper block 
in Fig. [H and simultaneously set in the range 8q G (— 90°, 0°) for the lower one. For the 
case 6q = 90° or 6q = 0°, this mode is identical to that of the single block in one unit cell, 
merged to the mode 2 or mode 4. The resonance modes predicted by analytical methods 
seems to agree quite well with the finite element simulation results. 

For different conjugate orientation angles, we obtained continuously changing resonance 
frequency under this coupling resonance mode (mode 3), with the other three resonance 
modes remains the same. Line with solid square in Fig. [5] shows that the resonance frequency 
dependant on the orientation angle of the elliptic cylinder. It decreases from 0°, and reaches 
the minimum at 90°. At both end of the curve, mode 3 disappeared, degenerating into 
that of the single block unit cell case. It also reveals that near 45°, the slope achieves the 
maximum value. Line with solid trigon in Fig. [5] shows the doubled block in unit cell case. 
With the same orientation angle 8q = 45° for both upper and lower blocks, the resonance 
frequency still dependent on the orientation angle, however, the vibration mode is no longer 
the same. This also lead to a different transmission property in the following section. 

B. Attenuation performances 

In order to provide efficient and reliable simulations for comparison with experimental 
data and other applications, the transmission properties is needed^ and calculated as fol- 
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FIG. 5: Orientation angle dependence of resonance frequency for conjugate blocks and double 
blocks in unit cell, respectively. 

lows. 

A periodic arrangement of unit cells with two coupling blocks illustrated in Fig. [1] is 
submerged in air. An incident plane wave is traveling along the positive x direction. Upper 
and lower boundaries were still exerted periodical boundary conditions. On left and right 
side, the air domains were added, governed by the equation 

V(-Vp) + — 2 v = 0, (20) 

which is provided by Time-harmonic analysis Module in COMSOL. The incident plane wave 
will encounter the air-epoxy boundary, inducing the displacement field in solids. The dis- 
placement field in solid domain is solved by employing Frequency response analysis Module. 
In order to illustrate the similarity of the equation form in the solid and air, we write the 
governing equations in the following form, 

where c, c p and c s are the wave speed in air, P-wave speed and 5*-wave speed in solid, 
respectively, and the repeated indices are summed over 1 to 3. These two modules are 
then coupled through domain boundaries by domain variables p and u. The harmonic 
acoustic pressure in the air on the interface acts as a boundary load pn to the solids. The 
model calculates harmonic displacements and stresses in the solids, and then it uses the 
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FIG. 6: Wave transmission spectrum in the low frequency region obtained from FEM simulations 
of a single layer composed by unit cells, the frequency resolution is 1Hz. Solid and dashed lines 
represent unit cell with conjugate blocks and double blocks, respectively. 

normal acceleration n • ii of the solid surface for air domain boundary to ensure continuity 
in acceleration. 

Fig. [6] shows the elastic wave transmission spectrum from 100Hz to 500Hz. Single 
frequency harmonic waves were employed in the simulation, and the minimum step is 1Hz. 
In the calculation, the absorption coefficient of the materials was taken to be zero, and 
transmission coefficient calculated for intensity of waves. The two dips and one peak observed 
in this figure are associated with the motion of inner resonator. The results from modal 
analysis are in good agreement with this transmission spectrum, in terms of the frequency 
positions of the two dips and the peak. The two dips, appearing at 313Hz and 363Hz, are 
corresponding to resonance vibration in Mode 2 and Mode 4. Mode 1 can not be observed in 
this figure, because the resonator remains symmetric in some extent. In the special case of 
elliptical resonators, whose center of mass away from its geometric center or the main axes 
are off the direction of wave propagation, rotational resonance can be expected to activate, 
bringing with an attenuation dip. For Mode 3, with conjugate orientation angles 45° and 
—45°, a peak at 332Hz is observed. This peak appears as an manifestation of the zero 
effective mass density for the block in x direction. The peak value is in fact 1, the lower 
value in the figure is an artifact due to the larger bin size in frequency used in the calculation. 
Accurate results around the resonance frequency will be given below. 

No phase shift can be observed in the transmission spectrum. The two dips in trans- 
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FIG. 7: Enlargement of the peak of the transmission spectrum in Fig. [6l the solid line represents 
the transmission peak in frequency step 0.01Hz, the circles are the fits to the calculated data with 



mission spectrum have been discussed by C. Goffaux— and M. Hirsekorn^. The resonance 
transmission peak and its corresponding mode is new. 

At the resonance frequency the peak value is close to 100%. Away from the resonance 
frequency the transmittance peak can be fitted neatly by the Lorentz form 



where uoq is the resonance frequency, F is the full width at half maximum of the resonance, 
and to the normalization constant. In Fig. [7] we plot the transmittance around the resonance 
frequency and the fit to a Lorentz curve, in this plot the frequency step is taken to be 0.01Hz. 
This result shows that the slab is effectively absent to the incident wave at a frequency near 
332Hz, as predicted at the end of Sec. [Til We denote this frequency by f z . It is worth noting 
that the f z is a little higher than f c , which in turn is slightly higher than the resonance 
frequency 332Hz. And all of the three are located in open interval (332Hz, 333Hz). When 
gradually tend to f z from the resonance frequency, p x e turns to zero, in other words, the 
effective thickness of the PC layer becomes smaller until zero, so that the attenuation effect 
is weakened. 

In comparison, we also calculated the transmittance of a system with a unit cell where 
the two elliptic cylinders in the same orientation (9 = 45°), which is a case of double blocks 
in unit cell of a simple PC. The result is shown in Fig. [6] as dashed line. It can be seen that 



Eq. {22]). 




(22) 
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the resonance modes 2 and 4 disappear completely, and mode 3 at frequency 332Hz (the 
cross point in Fig. [5]) turns to be a transmission dip. Resonance mode analysis indicates 
that the whole rigid matrix takes part in to partially balance the momentum in y direction 
induced by the cylinders, while cylinder resonance mainly exits in x direction. This lead 
to large EMD in x direction, and hence significant reflection for non-dissipative material. 
Actually, the minimum repeating unit of the system is the single cell. As investigated by 
Hirsekorn et a&, it makes no difference when the unit cell is doubled. 



IV. CONCLUSION 



In summary, the influence of the local asymmetry on the properties of the PCs at the low 
frequency in 2D has been studied. The analytical method used here are able to identify ro- 
tational mode, and the new mode whose vibration is perpendicular to the incident direction. 
We compared the results with the ones from finite element method, showing good agreement 
with each other. Transmission characteristic of a single layer PC are also calculated. The 
two dips are elastic wave attenuation corresponding to the common resonance modes, while 
the peak corresponds to the new resonance mode. Moreover, the peak can be tuned by 
simultaneously rotating the orientation angle of elliptical resonators of coupling blocks in a 
unit cell. This property provides a simple way to tune the resonance transmittance in locally 
resonant composite materials, leading to further opportunities for practical applications of 
the PCs. 



Electronic address: luoxd@sjtu.edu.cn 



1 Z. Liu, X. Zhang, Y. Mao, Y. Y. Zhu, Z. Yang, C. T. Chan, and P. Sheng, Science 289, 1734 
(2000). 

2 M. Hirsekorn, Appl. Phys. Lett. 84, 3364 (2004). 

3 M. Hirsekorn, P. P. Delsanto, A. C. Leung, P. Matic, J. Appl. Phys. 99, 124912 (2006). 

4 L. Y. Wu and L. W. Chen, J. Phys. D: Appl. Phys. 40, 7579 (2007). 

5 F. Cervera, L. Sanchis, et al, Phys. Rev. Lett. 88, 023902 (2002). 

6 C. Goffaux, F. Maseri, J. O. Vasseur, B. Djafari-Rouhani, and P. Lambin, Appl. Phys. Lett. 83, 
281 (2003). 

16 



J. V. Sanchez-Perez, C. Rubio, R. Martinez-Sala, R. Sanchez Grandia, and V. Gomez, Appl. 
Phys. Lett. 81, 5240 (2002). 

C. Goffaux and J. Sanchez-Dehesa, Phys. Rev. B 67, 144301 (2003). 

G. Wang, J. Wen, Y. Liu, and X. Wen, Phys. Rev. B 69, 184302 (2004). 

G. Wang, X. Wen, J. Wen, L. Shao, and Y. Liu, Phys. Rev. Lett. 93, 154302 (2004). 

For a review of this method see M.S. Kuswaha, Recent Res. Devel. Appl. Phys. 2, 743 (1999). 

M. Kafesaki, R. S. Penciu, and E. N. Economou, Phys. Rev. Lett. 84, 6050 (2000). 

M. G. Imhof, J. Acoust. Soc. Am. 100, 2969 (1996). 

Y. W. Gu, X. D. Luo, and H. R. Ma (to be published). 

M. G. Imhof, Geophys. J. Int. 156, 287 (2004). 

R. Egli, A. Geiger, A. Wiget and H. G. Kahle, Geophys. J. Int. 168, 1 (2007). 

A. Khelif, P. A. Deymier, B. Djafari-Rouhani, J. O. Vasseur and L. Dobrzynski, J. Appl. Phys. 

94, 1308, (2003). 

C. Goffaux, J. Sanchez-Dehesa, A. L. Yeyati, P. Lambin, A. Khelif, J. O. Vasseur, and B. 

Djafari-Rouhani, Phys. Rev. Lett. 88, 225502 (2002). 

Z. Liu, C. T. Chan, and P. Sheng, Phys. Rev. B 71, 014103 (2005). 

K.F.Graff, Wave Motion in Elastic Solids (Dover, New York, 1991). 

M. G. Imhof, J. Acoust. Soc. Am. 97, 754 (1995). 

A. J. Kalinowski, Except from the Proceedings of the COMSOL Conference 2007, Boston. 



17 



